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Abstract 

Measurements of the ZZ production cross sections in proton-proton collisions at 
center-of-mass energies of 7 and 8 TeV are presented. Candidate events for the 
leptonic decay mode ZZ —> 212v, where 1 denotes an electron or a muon, are 
reconstructed and selected from data corresponding to an integrated luminosity 
of 5.1 (19.6) fb^^ at 7 (8) TeV collected with the CMS experiment. The measured 
cross sections, cr{pp ZZ) = S.l+j'^ (stat) ±0.1(lumi)pb at 7 TeV, and 

7.2+q g (stat) + j'g (syst) ± 0.2 (lumi) pb at 8 TeV, are in good agreement with the stan¬ 
dard model predictions with next-to-leading-order accuracy. The selected data are 
analyzed to search for anomalous triple gauge couplings involving the ZZ final state. 
In the absence of any deviation from the standard model predictions, limits are set 
on the relevant parameters. These limits are then combined with the previously pub¬ 
lished CMS results for ZZ in 41 final states, yielding the most stringent constraints on 
the anomalous couplings. 
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1 Introduction 

The production of pairs of Z bosons in proton-proton collisions is a rare diboson process in 
the Standard Model (SM). The measurement of the cross section and properties of this process 
probe the self-interaction of electroweak gauge bosons. The ZZ final state is also an important 
background in searches for other interesting processes beyond the SM, such as the production 
of high-mass Higgs bosons and their subsequent decay to pairs of bosons fH or supersymme¬ 
try [[2j|. Because of the non-Abelian structure of the electroweak gauge theory, vector bosons 
can interact among themselves and can couple in triplets (e.g. WWZ) or quartets (e.g. WWZZ). 
All couplings involving only bosons without electric charge are expected to be null at tree level, 
leading to the absence of triple gauge couplings for Z 77 , ZZ 7 , and ZZZ. An enhancement in 
the measured rate of ZZ production compared to the expectation from the SM could indicate 
the existence of anomalous boson couplings. 

This paper presents measurements of the ZZ production cross sections in proton-proton col¬ 
lisions at the LHC at two different center-of-mass energies, 7 and 8 TeV, in the decay channel 
with two charged leptons, electrons (ee) or muons and a neutrino-antineutrino pair of 
any flavor (vv). The data were collected with the CMS detector at 7 ( 8 ) TeV, corresponding to 
5.1 (19.6) fb~^ of integrated luminosity. 

At tree level, ZZ pairs are primarily produced in the SM via the t- and w-channels, following the 
annihilation of a quark-antiquark pair in proton-proton collisions. Because of the high gluon- 
gluon parton luminosity, the gg —> ZZ contribution has to be included. The production cross 
section calculated up to next-to-leading-order (NLO) accuracy in strong coupling constant (as) 
is expected to be 6.46^^'^ (7.92+g 24 ) pb at 7 ( 8 ) TeV |3i|, where the uncertainties refer only to the 
missing higher orders in the computation. These cross sections include a leading-order (LO) 
computation of the gg —> ZZ contribution, which is formally a next-to-next-to-leading-order 
(NNLO) correction. Recently, complete NNLO cross sections for ZZ production accompanied 
by jets have also been computed |41, leading to a further small increase in cross section com¬ 
pared to Ref. |31. However, higher-order QCD corrections have been shown to be reduced 
significantly when vetoing events where the diboson system is produced in association with 
jets ||5l 13, as done in the present analysis. The NNLO QCD corrections apart from the LO 
gg — > ZZ contribution are thus neglected in our simulations and in the reference cross sections 
to which our measurements are compared. Complete one-loop electroweak (EW) corrections 
to massive vector boson pair production IZl | 8 l have also been published. The consequences of 
the EW corrections for ZZ production are that the transverse momentum (px) spectrum of the 
Z bosons falls more rapidly and, in addition, the overall cross section decreases by about 4% at 
EHC center-of-mass energies. 

The production of ZZ pairs has been studied at the EHC by the ATEAS experiment, which 
analyzed the decay modes 2121' and 212v (1,1' = e, p) at 7 TeV BH/ and by the CMS experiment, 
which considered 2121' final states (1 = e, p and 1' = e, p, t) at 7 TeV IflOl and 8 TeV IHtI . Both ex¬ 
periments measured ZZ production cross sections in good agreement with the SM predictions 
and set limits on anomalous triple gauge couplings (ATGCs). 

The branching fraction for the 212v decay mode (where 1 denotes only e and p) is approxi¬ 
mately six times larger than that of the four-charged-lepton final state, and the signal purity 
is enhanced at large values of the boson pj, where there is the greatest sensitivity to ATGC 
effects. Eor this reason, the 2l2v channel has a sensitivity comparable to that of the 41 channel 
to ATGG. The characteristic signature is an overall imbalance in the transverse momentum of 
the event between the initial and the final states, which consequently appears as missing trans¬ 
verse energy (£“^®®) in the final state. Although the branching fraction is large, this channel is 
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rather challenging due to the large contamination from background processes, in particular the 
Drell-Yan (DY) process, which has a cross section nearly five orders of magnitude larger than 
the signal. If the Z boson or the hadrons recoiling against it are not reconstructed correctly, then 
an apparent results and these events can resemble the signal. Other important sources of 
background are diboson processes, WW and WZ, with fully leptonic decays, and tf production. 

This paper presents a measurement of the ZZ production cross section in the 2l2v channel as a 
function of the transverse momentum (pj) of the charged lepton pair. The distribution of the 
dilepton pj is sensitive to the presence of ATGCs. Limits are computed and finally combined 
with existing results from CMS in the four-charged-lepton final state. 


2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam¬ 
eter, providing a magnetic field of 3.8 T. Within the superconducting solenoid volume are a 
silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and 
a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap 
sections. The silicon tracking system is used to measure the momentum of charged particles 
and covers the pseudorapidity range \rj\ < 2.5, where rj = —\n (tan (0/2)), and 0 is the po¬ 
lar angle of the trajectory of the particle with respect to the counterclockwise-beam direction. 
The ECAL and HCAL extend to a pseudorapidity range of |7/| < 3.0. A steel/quartz-fiber 
Cherenkov forward detector extends the calorimetric coverage to I?/! < 5.0. Muons are mea¬ 
sured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. 

The is defined as the magnitude of the missing transverse momentum or momentum 
imbalance, p^'®®, which is the negative vector sum of the momenta in the plane transverse to the 
beam of all reconstructed particles (photons, electrons, muons, charged and neutral hadrons) 
in the event. 

A more detailed description of the CMS detector, together with a definition of the coordinate 
system used and the relevant kinematic variables, can be found in Ref. IIT^ . 


3 Simulation 

Several Monte Carlo (MC) event generators are used to simulate the signal and background 
processes. The ZZ —)■ 2\2v signal and the WW —?■ 2\2v and WZ —> 31v background processes 
are simulated using MadCraph 5 fl3l . as well as Z+jets, W+jets, and tf+jets processes. Single 
top-quark processes are simulated with powheg IflTl . In the simulation, vector bosons are 
allowed to decay to leptons of any flavor (e, p, t), since t leptons can contribute to dielectron 
and dimuon final states through t — > e and t —> ^ decays. Eor all these processes, the parton 
showering is simulated with PYTHIA 6 1151 with the Z2 (Z2*) tune for 7 (8) TeV simulations 116] . 

The cross section of the ZZ signal is computed with the NLO generator mcfm |3l, which in¬ 
cludes contributions from gluon-gluon initial states. Since the present cross section measure¬ 
ment and ATCC analysis rely on the pj distribution of Z bosons, a precise prediction of this 
distribution is required. The charged dilepton pj spectrum of ZZ —> 212 p, generated with 
MadCraph and interfaced with PYTHIA for parton showering, is found to be in good agree¬ 
ment with the corresponding spectrum computed at NLO in QCD with MCFM and therefore no 
differential correction for NLO QCD effects is applied to the MadGraph simulated sample. 
In addition, the effect of NLO EW corrections (Zl El is taken into account by reweighting the 
ZZ and the WZ events as a function of the partonic kinematic variables, and applying weights 
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derived from the calculations described in Ref. [Zl. These corrections yield an overall reduction 
of 4.1% of the ZZ cross section, as well as a softening of the boson pj spectra that results in a 
reduction of the differential cross section of about 20% at Z pj of 300 GeV. 

Simulated samples of the ZZ —> 212v process that include contributions from ATGCs (see Sec¬ 
tion]^ are produced using the LO generator SHERPA HTtI . These samples are based on a LO 
matrix-element simulation including up to two additional jets, matched to parton showers. 

The parton distribution functions (PDF) are modeled with the GTEQ6L lllSl parametrization in 
samples generated with MadGraph, and the GTIO parametrization |[I9]| in samples generated 
with POWHEG and SHERPA. The detector response to the simulated events is modeled with 

geant 4 iniiza. 

4 Event selection 

The signal consists of two Z bosons, one decaying into a pair of oppositely charged leptons 
and the other to two neutrinos that escape direct detection. The final state is thus characterized 
by: a pair of oppositely charged, isolated electrons or muons, with an invariant mass within a 
Z-boson mass window, no additional leptons, and large 

Events are selected using triggers that require the presence of two electrons or two muons, with 
minimum pj thresholds on each lepton that depend on the dataset. The trigger thresholds in 
the 8TeV dataset are 17 and 8 GeV for the leptons with higher and lower pj, respectively. The 
thresholds for the 7 TeV data samples are the same or lower. The 8 TeV data sample also includes 
events that satisfy a single isolated muon trigger to ensure the highest efficiency. Eor events 
with two identified and isolated leptons having invariant mass between 83.5 and 98.5 GeV and 
dilepton pj > 45 GeV, the trigger efficiency is higher than 98% in the dielectron channel and 
varies from 94 to 98% in the dimuon channel. In addition, single-photon triggers or electron- 
muon triggers are used to select control samples for the background determinations. 

Electrons are selected inside the fiducial region of EGAE. The electron candidates must have 
a minimum pj of 20 GeV, and satisfy standard identification criteria, based on shower shape, 
track quality, cluster track matching, in order to reject misidentified hadrons |2TI. 

The muons are selected inside the fiducial region of the muon spectrometer, with a minimum 
Pj of 20 GeV, and satisfy standard identification criteria based on track information and isola¬ 
tion 11221. 

Events are selected if they include a pair of same-flavor, oppositely charged leptons that pass 
the identification and isolation criteria. In order to suppress backgrounds that do not include 
a Z boson, the lepton pair is required to have an invariant mass compatible with the Z-boson 
mass, between 83.5 and 98.5 GeV. The pj of the dilepton pair is required to be greater than 
45 GeV. This requirement is particularly effective at reducing the DY background because the 
Z bosons produced in ZZ events have, on average, larger pj than those from single Z-boson 
production. 

Since the ZZ pair is produced in the collision of two hadrons, the event might have jets from 
initial-state radiation. We use jets reconstructed from particle-flow (PE) candidates, using the 
anti-kx algorithm 1231 with a distance parameter of 0.5. The jet transverse energy is corrected 
using the GMS standard prescriptions for jet energy scale (JES) calibration Il24l . Only jets with 
a corrected pj greater than 10 GeV and reconstructed within |j/| <5 are used in this analy¬ 
sis. Eurther corrections are applied to reduce the effect of secondary proton-proton collisions 
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overlapping with the primary interaction (pileup). An extra correction is applied to jets in 
the MC samples to match the resolution observed in data. In order to reject jets dominated 
by instrumental and beam-related noise, loose identification criteria are applied, based on the 
multiplicity and energy fraction of charged and neutral particles. 

In order to suppress background coming from top quarks, events are vetoed if they have a jet 
identified as a b-quark jet (b-tagged). A requirement based on a combined secondary vertex 
discriminator Il25l is applied to b-tagged jets with pj > 20 GeV within the tracker fiducial 
region (|j/| < 2.4). The misidentification probability for light-parton jets is about 10%, whereas 
the efficiency for b-jets is more than 80%. To further reduce top-quark and other backgrounds 
with hadronic activity, events are rejected if they contain any jet with pj > 30 GeV. 

A good measurement is critical for the extraction of the ZZ —> 212v signal given that the 
pmiss distinguishes this process from the DY background. Since the average of the signal 
is moderate (~50GeV), we cannot simply require a high We follow the approach of 

constructing a "reduced variable, as done in the DO [2611271 and OPAL l(28l experiments. 

The concept behind a reduced is to reduce the instrumental contribution to mismeasured 
£miss considering possible contributions to fake In each event, and jet momenta 
are decomposed along an orthogonal set of axes in the transverse plane of the detector. One 
of the axes is defined by the px of the charged dilepton system, the other perpendicular to it. 
We define the recoil of the 1+D system in two different ways: (1) the clustered recoil (Rc) is the 
vectorial sum of the momenta of the PF jets reconstructed in the event, and (2) the unclustered 
recoil (Ru) is the vectorial sum of the transverse momenta of all PF candidates in the event, 
with the exception of the two leptons. On each axis (z = parallel/orthogonal to the dilepton 
system px), the reduced projection is defined as 

reduced =-p^Y - K/u^ 

where represents the choice of Rq or Ru that minimizes the absolute value of that reduced 

E^^®® component, and pY is a projection of the transverse momentum of the Z boson. The pres¬ 
ence of genuine E^^®® in the recoil of the charged dilepton system is expected to be evident in 
the parallel projection, while the component perpendicular to the 1+D system is mostly domi¬ 
nated by jet and E™^®® resolution. The absolute reduced E^^®® variable is the sum in quadrature 
of the two components. The reduced E“'®® shows better DY background suppression than the 
standard PF E^^®® at the same signal efficiency. It is also found to be more stable than the PF 
pmi®® variations in pileup conditions and JES. 

The E^^®® balance variable is defined as the ratio between the PF E3^®® and the transverse mo¬ 
mentum of the leptonically decaying Z boson, namely E^^®®/ pj. Values of this variable far from 
unity identify events in which the leptonic Z-boson candidate is not well balanced by genuine 
pmi®® neutrinos, but recoils against mismeasured jets or leptons. The selected sample can 
still be contaminated by events with jets with pj below the veto threshold. 

A mismeasurement of the jet energy can produce mismeasured p^^®® aligned with the jet di¬ 
rection in the transverse plane. These events are characterized by a small azimuthal angle 
between the p^®® vector and the closest jet, A^(p^®®, jet). This distribution is used to reject 
Z+jets events that have a small Acp angle. The mismeasurement of a lepton pj can also produce 
mismeasured E^^®®. Although this effect is usually negligible, given the good lepton momen¬ 
tum resolution in GMS, events are found where a large E^^®® value (>60GeV) is accompanied 
by a small angle between the p^'®® and the px of a lepton. Events with E“^®® > 60 GeV and 
A(/)(p^^®®, lepton) < 0.2 rad are therefore rejected. 
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In order to suppress the WZ background, with both bosons decaying leptonically, events are 
required to have no additional leptons. To improve the rejection power, the px threshold is low¬ 
ered to 3 GeV for additional muons, and 10 GeV for electrons. Furthermore, these muons and 
electrons are selected with looser criteria than those used to reconstruct the Z-boson candidate. 

The variables described above are used to extract the signal sample for the cross section mea¬ 
surement. We optimize the requirements in the final selection in order to minimize the total 
uncertainty in the measured cross section at 8 TeV (see Section]^. The same selection is applied 
to the 7 TeV data. For this purpose, we scan a series of possible analysis selections, in which 
we vary the dilepton mass window and px threshold, the minimum px of jets used in the com¬ 
putation of the reduced variable, and the reduced requirement. We optimize the 
selection using MG estimates of the background processes, or using predictions based on con¬ 
trol samples in data from the DY, top-quark, and WW backgrounds, as described in Section]^ 
and we find similar results for the optimal requirements and for the measured cross section. 
For the final optimization we choose the selection obtained using background estimates from 
data. The requirements are summarized in Table With this selection, the acceptance for 
ZZ —> 2e2v and ZZ —)■ 2p2t/ events is about 10% for both channels, at 7 and 8 TeV. 


Table 1: Summary of the optimal signal selection. 


Variable 

Value 

Dilepton invariant mass 

|m(ll) -91| < 7.5 GeV 

Dilepton pj 

pl^ > 45 GeV 

b-tagged jets 

Based on vertex info (for jet with pj > 20 GeV) 

Jet veto 

No jets with px > 30 GeV 

Reduced £J^'®® 

>65 GeV 

£miss tjalance 

0.4 < £^i®®/p!jl < 1.8 

A^(pj‘®®, jet) 

>0.5 rad 

A^(p!JT®®, lepton) 

> 0.2 rad 

Lepton veto 

No additional leptons (e/p) with pj > 10/3 GeV 


5 Background estimation 

Although the DY process does not include genuine from neutrinos, the tail of the reduced 

distribution can be contaminated by these events due to detector energy resolution, jet 
energy mismeasurements, pileup energy fluctuations, and instrumental noise. Given that the 
simulation may not fully reproduce detector and pileup effects on the reduced £^^®® distribu¬ 
tion, especially in the tails, and that the simulation is limited in statistical precision, we build a 
model of DY background from control samples in data. For this purpose we use a process that 
has similar jet multiplicity, underlying event, and pileup conditions as the DY process for the 
region of interest at high boson pj: the production of prompt isolated photons in association 
with jets (7 + jets) Il29l . We expect that an accurate description of the £^^®® distribution and 
other related kinematic variables can be obtained from this photon + jets sample. However, 
some corrections must be applied to the photon + jets sample to ensure a good modeling of the 
DY process. The yield of photon events is scaled to the observed charged dilepton system yield 
as a function of the boson pj after applying the jet veto to both samples. This accounts for the 
differences in the selection efficiency of the dilepton and photon candidates and corrects for the 
trigger prescales, which are applied to the low-px photon triggers. 
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Only photons in the barrel region are used because the purity and resolution are better than 
in other regions. Following Ref. [1], the selection of photon events is based on shower shape, 
isolation in the tracker, and energy deposits in ECAL, and HCAL. After this selection, several 
processes with instrumental contribute to the photon sample: single 7 events, double 7 
events where one photon escapes detection or fails the identification, and QCD events with a 
jet misidentified as a photon. Processes with genuine can also contaminate this sample: 
W/Z +7 with the W/Z boson decaying to \vIvv, or W+jets with the W boson decaying to ev 
and the electron misreconstructed as a photon. Although these processes have generally lower 
cross sections, they are characterized by large values, and thus contribute to the tails of 
the distribution, where it is most important to measure the residual instrumental background. 
In order to reduce these background contributions, specific selections are applied. The event 

must have exactly one photon and no leptons. Only jets with AR = \/{A(py - + {Ailf > 0.4 
from the photon are used for all the jet-related selections (jet veto, reduced £^'®®, etc.). To 
avoid misreconstruction of the photon energy, a conversion veto is applied using the number 
of missing expected tracker hits and the distance of closest approach between the reconstructed 
conversion tracks. 

The remaining contribution from W + 7 and W/Z -|- 7 events after this selection is estimated 
from simulation and subtracted from the photon data model. For this purpose, a set of sim¬ 
ulated photon samples is used that includes 7 + jets, QCD events with a jet misidentified as 
a photon (generated with pythia), W -|- 7 —> \vj, and Z -|- 7 —)■ 1 / 1/7 (generated with Mad- 
Graph). These samples are normalized to their respective cross sections computed at NLO in 
QCD. The full set of MC samples is reweighted and corrected following the same procedure as 
that used for the photon data sample. Finally, the photon data are corrected as a function of 
£miss multiplying them by unity minus the fraction of electroweak processes in the simula¬ 
tion. 

We apply a different data-based method to estimate the total number of background events 
from processes that do not involve a Z boson: i.e. WW and top-quark production. We denote 
these events as nonresonant background (NRB). In order to measure this contribution, a control 
sample based on e/i candidate events is selected by applying the same requirements as in the 
main analysis. The NRB yields in the same-flavor charmels (ee and jiji) are obtained by scaling 
the number of events in the control sample. The rescaling is done by means of correction 
factors, measured from the sidebands (SB) of the Z-boson mass peak, i.e. in the regions 55- 
70 and 110-200 GeV. The scale factors are measured in a looser selection region in order to 
improve the statistical precision. We require the reduced £^^®® > 65 GeV in order to suppress 
the DY contribution from t+t~. We also require at least one b-tagged jet with pj > 20 GeV, to 
further reduce DY and other backgrounds, and increase the fraction of top-quark events. The 
scale factors are defined as follows: 


^ _ XfSB / ArSB 


and the NRB contamination in the Z-peak region is: 




peak 
efi ■ 


( 1 ) 

( 2 ) 


The validity of the method is tested in simulation by comparing the predicted background to 
the expected number of WW and top-quark events. 

Figure [1 shows the reduced £“'®® distributions in dilepton data and simulation, using the pho¬ 
ton model to describe the DY background and the data-driven estimation for NRB. A good 
agreement is found in the region dominated by the DY process, up to about 80 GeV, while the 
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5.1 fb ' (7 TeV) 



Reduced E™"" [GeV] 


19.6fb ' (8 TeV) 



Reduced E™"" [GeV] 


Figure 1: Reduced spectrum in the inclusive 11 (1 = e, channel at 7 TeV (left) and 8 TeV 
(right), using the photon model to describe the DY contribution and NRB modeling for WW, 
W + jets, and top-quark production, after selections on the dilepton invariant mass and pj, jet 
veto, b-tagged jet veto, third lepton veto, and jet), as described in Section]^ The gray 

error band represents the statistical uncertainty in the predicted yields. 


higher part of the spectrum is dominated by diboson production. The error bands shown in 
Fig. represent the statistical uncertainty in the predicted yields. A systematic uncertainty 
in the final DY event yield estimated with this method is computed as the relative difference 
between dilepton yields in data and simulation, in a control region with < 60 GeV, and 
it has been found to be 25% (40%) at 7 ( 8 ) TeV. This systematic uncertainty is not shown in 
Figure 


6 Systematic uncertainties 

Different sources of systematic uncertainty are associated with the expected yields and distri¬ 
butions of signal and background processes and of the data. The uncertainties reported in the 
following paragraphs affect the final event yields of the relevant processes. 

Statistical uncertainty of the simulated and control samples. For the processes esti¬ 
mated from simulation, ZZ and WZ, the limited size of the MC sample affects the precision 
of the modeling, and is therefore taken as a systematic uncertainty in the shape of the kine¬ 
matic distributions used in the cross section measurement and ATGC limit setting. Similarly, 
the backgrounds estimated from data are limited by the size of the control samples described 
in Section]^ the ep sample for nonresonant backgrounds and the 7 +jets sample for DY back¬ 
ground. These uncertainties are treated in the same way as those backgrounds that are esti¬ 
mated from simulation. This systematic uncertainty has been computed in different reduced 
£miss Qj. cJifferent pj bins and is used as shape errors in the fit. 

Cross sections of ZZ and WZ. The cross sections for pp —> ZZ -|- X —> 212v + X and 
pp — > WZ -|- X — > 3lv + X processes are calculated using mcfm version 6.2 |3|, and using 
PDFs from the Les Houches accord PDF (lhapdf) program, version 5.8.7 K301I . The PDF+as 
uncertainty in the WZ cross section is evaluated as the maximum spread of the cross sections 
computed at Pr = pp = with three PDF sets, including the corresponding uncertainties 
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6 Systematic uncertainties 


from one standard deviation variation of the PDF parameters and the value (31]. It is found 
to be 3.1% (4.2%) at 7 (8) TeV. 

The uncertainty from the renormalization and factorization scales is evaluated as the maximum 
difference between the central value of the cross section at = jip = niz and the central val¬ 
ues computed at }Ir = fip = ?r?z/2 and 2niz, using each of the three PDFs recommended in 
Ref. IlSn . An uncertainty of 5.9% (5.4%) at 7 (8) TeV is found for the WZ background. For the 
ZZ signal, we evaluate this theoretical uncertainty in the case of the exclusive production with 
0 jets, to take into account the jet-veto applied in the signal selection, following the prescrip¬ 
tion described in Refs. ||32H^. The exclusive cross section for ZZ+0 jets is (Toy = — cr^y, 

where cry„j is the inclusive cross section of ZZ -|- at least n jets, where n = 0,1. According to 
Ref. Il32l , cr^oy and (7>iy are essentially uncorrelated, thus the uncertainty in tjQj can be computed 

as €oj = 'J e>oy + ^>iy/ where e>oy and 6>iy are the uncertainties in (r^oy and respectively. 
The cross sections are computed with MCFM, including the acceptance requirements on lepton 
Pj and rj, charged dilepton mass, and as well as the jet veto, when relevant. The cross 

section uncertainties are estimated by varying the renormalization and factorization scales, as 
explained above. Since the charged dilepton pj spectrum is the observable from which lim¬ 
its on ATGCs are derived, the uncertainty in (Joy is computed in different intervals of charged 
dilepton pj. 

The uncertainty in the NLO EW correction to ZZ production, corresponding to missing higher- 
order terms in the computation, is estimated as the product of the NLO QCD and EW cor¬ 
rections |7|. The uncertainty in the EW correction to WZ production is estimated as 100% of 
the correction, to account for the poorly known fraction of photon+quark-induced events Il8ll 
passing the jet veto. 

Acceptance. The kinematic acceptance for the signal is computed using MCFM. Kinematic 
requirements, based on those used in the signal selection, are applied to the charged leptons 
and neutrinos at the generator level. The acceptance is determined by comparing the cross sec¬ 
tions with and without the kinematic requirements. The systematic uncertainty is evaluated 
as the variation in the acceptance resulting from varying the renormalization and factorization 
scales from niz to mzl2 and 2 mz, summed in quadrature with the variation obtained from us¬ 
ing different FDE sets and from varying the FDE parameters and the value by one standard 
deviation. The result is 2.8% at both 7 and 8 TeV. 

Luminosity. The uncertainty in the luminosity measurement is 2.2% in 2011, and 2.6% in 
2012 1341 . 

Lepton trigger and identification efficiency. Lepton trigger and identification efficiencies 
are determined from data, using the tag-and-probe technique with Z —> 11 events Il35l . and 
used to correct the simulated samples. The total uncertainty in the lepton efficiency amounts 
to about 3% for ee events, and 4% for pp events. 

Lepton momentum scaie. The systematic uncertainty in the lepton momentum scale is 
computed by shifting the nominal momenta by ±1(J and propagating the variations to the 
reduced We assume an uncertainty of 2% (3.5%) in the energy of electrons reconstructed 
in the ECAL barrel (endcap), and 1% in the muon momentum. The resulting variations of the 
final yields are 2.5% for the ee channel, and 1.0% for the pp channel and they are treated as a 
shape uncertainty. 
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Jet energy scale and resolution. The uncertainty in the calibration of the jet energy scale 
directly affects the jet veto, the calculation of reduced and the selection of the balance 

variable. The JES uncertainty is estimated by shifting the jet energies by iltr and propagating 
the variations to the reduced and all the other relevant observables. Uncertainties in the 
final yields of 3-4 (7-8)% are found for both the ee and final states at 7 (8) TeV. 

Similarly, a systematic uncertainty in jet energy resolution (JER) is computed. As explained 
above, the energy of jets in simulation is corrected to reproduce the resolution observed in data. 
Such corrections are varied according to their uncertainties and these variations are propagated 
to all the observables and selections dependent on jet energy. An uncertainty in the final yields 
of less than 1% is found in both ee and jiji final states: 0.4% (0.8%) at 7 (8) TeV. 

Since the shapes of the distributions are expected to be affected by variations in the JES and the 
JER, these sources are treated as shape uncertainties in the extraction of the cross section. 

b-jet veto. The b-tagging efficiency is taken from Ref. I36l . In simulation, the nominal 
working point for this b-tagger is shifted to reproduce the efficiency observed in data. The 
uncertainty in the measured efficiency is propagated to the event yields of the processes esti¬ 
mated from simulation by applying further shifts to the discriminator threshold. A very small 
uncertainty in the final yields of the MC samples is found: 0.1-0.15% at both 7 and 8 TeV. 

Plleup. Simulated samples are reweighted to reproduce the pileup conditions observed in 
data. To compute the uncertainty related to this procedure, we shift the number of interactions 
by 8% when reweightmg the simulated samples. The variation of the final yields induced by 
this procedure is less than 1% in ZZ and WZ processes. However, the shapes of the kinematic 
distributions can vary in this procedure, so the varied distributions are used as shape uncer¬ 
tainties in the cross section fit. 

Drell-Yan. The uncertainty in the DY contribution is propagated from the uncertainty in 
the reweighted photon spectrum that is used in the estimate of DY background from data, and 
is dominated by the subtraction of backgrounds due to EW processes. As explained in Sec¬ 
tion 1^ the DY background estimate is assigned an uncertainty of 25% (40%) at 7 (8) TeV, eval¬ 
uated from the relative difference between dilepton yields in data and simulation in a control 
region. 

Top-quark and WW backgrounds. The uncertainty in the estimate of the NRB is derived 
from the statistical uncertainties in the scale factors in Eq. Q, and from a closure test of the 
data-driven method for the measurement of this background performed on simulated data. It 
is found to be about 20% at both 7 and 8 TeV. 


7 Measurement of the ZZ production cross section 

We extract the ZZ production cross section using a profile likelihood fit Il37l to the reduced-£^®® 
distribution, shown in Eig. The fit takes into account the expectations for the different back¬ 
ground processes and the ZZ signal. Each systematic uncertainty is introduced in the fit as a 
nuisance parameter with a log-normal prior. Eor the signal we consider a further multiplicative 
factor, which is the ratio of the cross section measured in data to the expected theoretical value, 
i.e. the signal strength }i = a/ c^th. Maximizing the profile likelihood, we obtain the ZZ produc¬ 
tion cross section from the signal strength parameter, as well as optimal fits of the background 
yields by varying nuisance parameters within their constraints. Table shows the expected 



10 


7 Measurement of the ZZ production cross section 


signal and background yields, and the corresponding values after the combined fit to the ee 
and fiji channels. The uncertainties include both the statistical and systematic components. 

Table 2: Predicted signal and background yields at 7 and 8 TeV, and corresponding values ob¬ 
tained from the combined maximum likelihood fit to the ee and channels. The uncertainties 
include both the statistical and systematic components. 


Dataset 

Process 

Channel 

Predicted yield 

Fitted yield 

Observed 


ZZ ^ 212v 

ee 

14.0 ± 1.9 

12.0 ± 4.4 

— 


m 

21.7 ±3.2 

18.4 ±6.8 

— 


WZ ^ 31v 

ee 

7.7 ±0.9 

7.9 ± 1.0 

— 


m 

11.5 ±1.6 

11.6 ± 1.2 

— 

7 TeV 

Z + jets 

ee 

5.0 ±2.7 

4.8 ± 2.3 

_ 

m 

8.3 ±4.8 

4.8 ± 3.0 

— 


Nonresonant 

ee 

7.7 ±3.1 

7.4 ± 2.3 

— 


W 

11.2 ±4.8 

9.2 ±3.1 

— 


Total 

ee 

34.4 ± 6.2 

32.1 ± 3.9 

35 


W 

52.7 ±9.7 

44.0 ± 5.3 

40 


ZZ ^ 212v 

ee 

77 ±16 

69 ±13 

— 


m 

109 ± 23 

100 ± 19 

— 


WZ ^ 31v 

ee 

45 ±6 

43.9 ± 5.6 

— 


m 

64 ± 8 

63.8 ± 7.3 

— 

8 TeV 

Z + jets 

ee 

36 ±12 

27.7 ± 7.9 

_ 

m 

63 ±21 

52 ±14 

— 


Nonresonant 

ee 

31 ±9 

34.1 ± 7.2 

— 


m 

50 ± 14 

54± 12 

— 


Total 

ee 

189 ± 31 

174.7 ± 10 

176 


W 

286 ± 49 

269.8 ± 15 

271 
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Figure 2: Reduced distribution in 11 (1 = e,}i) channels, after the full selection, at 7 TeV 
(left) and 8 TeV (right). The DY and WW, W+jets, and top backgrounds are estimated with 
data-driven methods. The gray error band includes statistical and systematic uncertainties in 
the predicted yields. In the bottom plots, vertical error bars and bands are relative to the total 
predicted yields. In all plots, horizontal error bars indicate the bin width. 
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The cross sections are extracted from individual fits to the ee and jiji channels and from a 
simultaneous fit to both channels. Table reports the measured pp —> ZZ —> 212v exclusive 
cross section, i.e. the production cross section of ZZ pairs with mass 60 < Mz < 120 GeV, with 
no restrictions on lepton acceptance nor jet number, times the branching fraction to final states 
with two charged leptons of a given flavor and two neutrinos of any flavor. This is obtained by 
rescaling the theoretical prediction for the exclusive cross section in the same kinematic range 
by the fitted signal strength. These theoretical predictions are computed at NLO in QCD with 
MCFM and corrected for NLO EW effects: 79^3 i97^f} fb at 7 (8) TeV. 

Table 3: Cross sections (fb) for process pp ZZ —> 212v (where 1 denotes a charged lepton of a 
given flavor, v a neutrino of any flavor) at 7 and 8 TeV, with both Z boson masses in the range 
60-120 GeV, measured in the ee and ]i}i channels and the two channels combined. 


Ghannel 

Vs = 7TeV 

Vs = 8 TeV 

ee 

pp 

98 I 33 (stat) I 22 ^ 2 (lumi) 

47 I 21 (stat) (syst) ± 1 (lumi) 

83l|g (stat) (syst) ± 2 (lumi) 

98 II 4 (stat) ^22 (sy®t) ± 3 (lumi) 

Gombined 

66^38 (stat) I]® (syst) ± 1 (lumi) 

92 ZIJ (stat) (syst) ± 2 (lumi) 

Theory 

79tt (theo) 

97tt (theo) 


The measured inclusive ZZ cross section is obtained by rescaling the theoretical inclusive cross 
section computed in the zero-width approximation |j3l and corrected for NLO EW effects (Zl 
(see Section [^, by the same fitted signal strength. This procedure properly accounts for the 
contribution of virtual photon decays to the charged-lepton pair production, and yields a mea¬ 
sured cross section that can be compared directly with theoretical calculations of inclusive pure 
ZZ production in the zero-width approximation. The results are: 

7TeV : cr(pp —> ZZ) = S.l^J'^ (stat) (syst) ± 0.1 (lumi)pb, 

8TeV : cr[pp —> ZZ) = 7.2^08 (stat) (syst) ± 0.2 (lumi)pb. 

This is the first cross section measurement in the 2\2v channel at 8 TeV. The measurements are 
less than one standard deviation from the SM predictions at both 7 and 8 TeV. The uncertainties 
are approximately twice as large as those from the GMS measurement in the 4£ channel nniiTii, 
and the charmels agree within uncertainties. 

The p-values of the simultaneous fit to the ee and pp channels are 0.335 (0.569) at 7 (8) TeV. The 
data are also consistent with the reduced spectra uncorrected for NLO EW effects, but 
with slightly smaller p-values of 0.322 (0.477) at 7 (8) TeV. The application of EW corrections 
thus improves the modeling of the diboson processes and leads to a better agreement between 
the simulated and observed spectra. 

Table 1^ shows a summary of the sources of systematic uncertainty described in Section]^ with 
the corresponding contributions to the total systematic uncertainty in the cross sections. 

8 Anomalous couplings 

The existence of neutral trilinear gauge couplings is forbidden at the tree level, but allowed 
in some extensions of the SM Il38l . The ZZ production process provides a way to probe the 
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8 Anomalous couplings 


Table 4: Systematic uncertainties in the cross sections due to each source separately, after the 
maximum likelihood fit to extract the ZZ cross section. The uncertainties marked with an 
asterisk (*) are used as shape uncertainties in the fit. 


Source of uncertainty 

Uncertainty [%] 
7TeV 8 TeV 

(*) MC statistics: ZZ (ee channel) 

0.8 

1.0 

(*) MC statistics: ZZ {jiji channel) 

1.3 

1.1 

(*) MC statistics: WZ (ee channel) 

1.7 

0.9 

(*) MC statistics: WZ {jiji channel) 

1.7 

1.0 

(*) Control sample statistics: DY (ee channel) 

6.9 

2.3 

(*) Control sample statistics: DY channel) 

5.8 

4.9 

(*) Control sample statistics: NRB (ee channel) 

6.3 

3.0 

(*) Control sample statistics: NRB channel) 

8.1 

4.4 

WZ cross section: PDP+ag 

1.9 

2.6 

(*) ZZ + WZ cross section: scales 

17 

16 

(*) ZZ + WZ cross section: NLO EW corn 

2.4 

2.3 

Signal acceptance 

2.8 

2.8 

(*) Pileup 

0.5 

1.1 

Muon trigger, ID, isolation 

4.1 

3.6 

Electron trigger, ID, isolation 

1.7 

2.0 

(*) Lepton momentum scale 

2.6 

3.7 

(*)JES 

6.0 

12 

{*) JER 

0.8 

1.4 

(*) Unclustered 

2.1 

3.2 

(*) b-jet veto 

0.3 

0.5 

Drell-Yan bkg. normalization 

6.6 

8.4 

Top-quark & WW bkg. normalization 

7.7 

7.1 

Total systematic uncertainty 

24.6 

23.5 

Statistical uncertainty 

28.0 

11.9 


existence of such anomalous couplings at the ZZZ and yZZ vertices. 

Neutral couplings V*'*^ZZ (V = Z, 7 ) can be described using the following effective Lagran- 
gian [391: 


-^vzz = “Tvy 
M| 


{[fl +/f - [f^ +/ 5 Z (a^^z,,)] 


(3) 


where Z represents the Z boson and Ff^ represents the electromagnetic field tensor. The co¬ 
efficients and /? correspond to couplings j^*'>ZZ and Z^*'>ZZ, respectively. All the oper¬ 
ators in Eq. 0 are Lorentz-invariant and U(1)em gauge-invariant, but not invariant under 
SU(2 )l X U(1)y gauge symmetry. The terms corresponding to fj parameters violate the CP 
symmetry, while the terms corresponding to /J parameters conserve CP. 


To avoid unitarity violation at energies above the scale (A) of new physics, the Lagrangian of 
Eq. (j^ can be modified with form factors of the type 1/ {1 + s/A)", where \/s is the effective 
center-of-mass energy of the collision. No form-factor scaling is used in this analysis. This 
allows to provide results without any bias that can arise due to a particular choice of the form- 
factor energy dependence. 

Previous studies of neutral anomalous triple gauge couplings were performed at LEP2 Il40l . 
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Tevatron HSl , and LHC EHTTI. No deviation from the SM expectation has been observed so 
far, and the best limits were set by the LHC measurements based on integrated luminosities of 
about 5 (19.6) fb^^ at 7 (8) TeV. 


8.1 Limits from the ZZ 2l2v channei 

In the following, we extract limits on the neutral triple gauge couplings V(*^ZZ with the same 
datasets at 7 and 8 TeV as used for the ZZ cross section measurement described in the previ¬ 
ous section. Limits on the four parameters are set by comparing the data with theoretical 
predictions. 

Figure|^shows the charged dilepton pj distribution after the full selection described in Table[^ 
in data and simulation, including SHERPA samples with different values of the parameter. 
The contribution from the anomalous couplings enhances the high-pj region of the distribu¬ 
tion. The charged dilepton pj is thus a good observable to probe for the presence of ATGCs. 
The DY and nonresonant backgrounds are estimated from data as described above. The SM 
ZZ process is simulated here using the MadGraph sample described in Section]^ with NLO 
QGD corrections computed with MCFM and NLO EW corrections from Ref. |ZJ- The contri¬ 
bution of the ATGGs is obtained from the SHERPA samples mentioned above, by subtracting 
the SM SHERPA contribution to the charged dilepton pj, and is summed to the MadGraph 
ZZ distribution. The interference of the ATGG signal and the SM ZZ production is included, 
except for Pt{Z) < 200 GeV, which has a negligible impact on the limits. The expected signal 
yields in each pj bin are interpolated between different values of the ATGG coupling parame¬ 
ters using a second-degree polynomial, since the signal cross section depends quadratically on 
such parameters. 
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Figure 3: Dilepton (1 = e, p) transverse momentum distributions at 7TeV (left) and 8 TeV (right). 
The DY and WW, W+jets, and top backgrounds are estimated from control samples in data. The 
gray error band includes statistical and systematic uncertainties in the predicted yields. In the 
bottom plots, vertical error bars and bands are relative to the total predicted yields. In all plots, 
horizontal error bars indicate the bin width. 


The limits are calculated with a profile likelihood method. We set one-dimensional limits on 
the four parameters, i.e. varying independently a single parameter at a time, while fixing the 
other three to zero. The 95% GL one-dimensional limits on the four parameters are reported 
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9 Summary 


in Table for 7 TeV, 8 TeV, and combined datasets. The observed exclusion limits are about 
one standard deviation tighter than the expected ones, which is attributed primarily to the 
observed deficit of events in the highest bin of dilepton px- The limits set are of comparable 
sensitivity to those previously obtained by CMS in the 41 channel IITOlITTl . 


Table 5: Summary of 95% CL intervals for the neutral ATGC coefficients, set by the IXlv final 
states using the 7 and 8 TeV CMS datasets. The expected 95% CL intervals obtained using the 7 
and 8 TeV simulated samples are also shown. No form factor is used. 


Dataset 

fi fl fl 

7 TeV 

[-0.010; 0.011] 

[-0.012; 0.013] 

[-0.010; 0.010] 

[-0.013; 0.013] 

8 TeV 

[-0.0033; 0.0037] 

[-0.0044; 0.0038] 

[-0.0033; 0.0035] 

[-0.0039; 0.0043] 

Combined 

[-0.0028; 0.0032] 

[-0.0037; 0.0033] 

[-0.0029; 0.0031] 

[-0.0033; 0.0037] 

Expected (7 and 8 TeV) 

[-0.0048; 0.0051] 

[-0.0060; 0.0053] 

[-0.0048; 0.0050] 

[-0.0057; 0.0062] 


8.2 Combined limits from the ZZ 41 and 2l2i/ channels 

We proceed with the combination of the results of the previously published ZZ —> 41 analy¬ 
ses Col HD with the present results. In doing this, the published analysis of the 41 (1 = e, p) 
channel is unchanged, except that NLO EW corrections to the SM ZZ —> 41 background are 
included in the same way as in the present analysis. We use a profile likelihood method to 
calculate the 95% CL one-dimensional intervals for the four parameters, combining the data in 
the 41 and IXlv channels, at 7 and 8 TeV. The systematic uncertainties in the signal and dibo¬ 
son background cross sections, in the integrated luminosity, and in the lepton efficiencies are 
treated as fully correlated between the two channels. Table shows the intervals obtained by 
combining the four separate data sets. The combined analysis improves the sensitivity of the 
two separate channels, and the limits are more stringent than all the results published to date. 

Table 6: Summary of 95% CL intervals for the neutral ATGC coefficients, set by the combined 
analysis of 41 and IXlv final states. The intervals obtained separately by the two analyses using 
the 7 and 8 TeV CMS data sets are shown, as well as their combination. The expected 95% CL 
intervals obtained using the 7 and 8 TeV simulated samples of both analyses are also shown. 
No form factor is used. 


Dataset 

fl 

fl 

fi 

r. 

7 TeV, 41 
7TeV,212v 

8 TeV, 41 
8TeV,212v 

[-0.010; 0.011] 
[-0.010; 0.011] 
[-0.0041; 0.0044] 
[-0.0033; 0.0037] 

[-0.012; 0.013] 
[-0.012; 0.013] 
[-0.0052; 0.0048] 
[-0.0044; 0.0038] 

[-0.011; 0.011] 
[-0.010; 0.010] 
[-0.0041; 0.0040] 
[-0.0033; 0.0035] 

[-0.013; 0.013] 
[-0.013; 0.013] 
[-0.0048; 0.0045] 
[-0.0039; 0.0043] 

Combined 

[-0.0022; 0.0026] 

[-0.0029; 0.0026] 

[-0.0023; 0.0023] 

[-0.0026; 0.0027] 

Expected 

(41 and 212v, 7 and 8 TeV) 

[-0.0036; 0.0039] 

[-0.0046; 0.0041] 

[-0.0036; 0.0037] 

[-0.0043; 0.0043] 


9 Summary 

We have measured the ZZ production cross section in the 212v channel in proton-proton colli¬ 
sions at center-of-mass energies of 7 and 8 TeV. The data samples selected for the study corre¬ 
spond to an integrated luminosity of 5.1 (19.6) fb~^ at 7 (8) TeV. We have measured 

c(pp —> ZZ) = 5.11J'4 (stat) (syst) ± 0.1 (lumi) pb 
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at 7 TeV, and 

cr(pp —> ZZ) = 7.2+Qg (stat) (syst) ± 0.2(lumi)pb 

at 8TeV, in agreement with theory calculations, 6.2^q\ pb (Z.blggpb) at 7(8) TeV, which include 
NLO QCD corrections iQ and NLO EW corrections IZlIBl- The selected data have also been 
analyzed to search for ATGCs involving the ZZ final state. In the absence of any observation of 
new physics, we have set the most stringent limits to date on the relevant ATGC parameters. In 
addition, by combining the selected data with the GMS data for the four-charged-lepton final 
state we have set even tighter constraints. 
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